Anharmonic vs. relaxational sound damping in glasses: 

II. Vitreous silica 



in 
o 
o 

(N 
> 

o 

(N 
(N 

C 

CZ5 



C3 



I 

C 

o 
o 



(N 
> 

o 

in 
in 
o 
in 
o 



c 

o 
o 



X 
S3 



Rene Vacher, Eric Courtens, and Marie Foret 
Groupe de Physique des Verres et Spectroscopies, LCVN, UMR CNRS 5587, 
Universite Montpelher II, F-34095 Montpellicr Cedex 5, France 
(Dated: February 2, 2008) 

The temperature dependence of the frequency dispersion in the sound velocity and damping of 
vitreous silica is reanalyzed. Thermally activated relaxation accounts for the sound attenuation 
observed above 10 K at sonic and ultrasonic frequencies. Its extrapolation to the hypersonic regime 
reveals that the anharmonic coupling to the thermal bath becomes important in Brillouin-scattering 
measurements. At 35 GHz and room temperature, the damping due to this anharmonicity is found 
to be nearly twice that produced by thermally activated relaxation. The analysis also reveals a 
sizeable velocity increase with temperature which is not related with sound dispersion. A possible 
explanation is that silica experiences a gradual structural change that already starts well below 
room temperature. 

PACS numbers: 63.50.+X, 78.40.Pg, 62.80.+f, 78.35,+c 



I. INTRODUCTION 

In the first paper of this series (I), the hypersonic atten- 
uation of highly dcnsificd silica glass, <i-Si02, was investi- 
gated. It was found that in that material the damping of 
hypersound is completely dominated by its anharmonic 
interaction with the thermally dominant modes. In the 
present paper we consider the corresponding situation in 
usual vitreous silica, v-Si02 , for which a large quantity of 
high quality data is already available in the literature. In 
the early 1950's, ultrasonic absorption peaks in function 
of the temperature T were first observed in z;-Si02- These 
were described by Anderson and Bommcl in terms of a 
phenomenological Maxwell model pj. Following the dis- 
covery of the key role played by two-level systems (TLS) 
in producing the anomalous thermal properties of glasses 
it became clear that these should also be invoked 
in the description of acoustic relaxation. A theory in- 
cluding both resonant and relaxational sound damping 
produced by tunneling was developed by Jackie [^|. It 
was then extended to higher temperatures by including 
in the description the classical jumps over the energy 
barrier separating the wells of TLS [3 . The presently ac- 
cepted model for thermally activated relaxation (TAR) 
includes a distribution of asymmetric double-well poten- 
tials, as discussed by Gilroy and Phillips |fj. 

The simple extrapolation of this behavior to GHz fre- 
quencies, i.e. to the hypersonic regime, has sometimes 
been presented as an appropriate description of the ob- 
served phenomena, e.g. in 0,0] or Q. However, it has 
also been recognized that TAR might not always be suf- 
ficient to account for sound dispersion and damping. A 
specific example has been discussed in the case of vitre- 
ous germania, GeC>2 The results presented in (I) 
strongly suggest that v-Si02 should be an excellent can- 
didate to search for a possible anharmonic contribution 
to the hypersonic attenuation. 

The paper is organized as follows. In Section II the 
phenomenology of TAR is reviewed to formulate a de- 



scription that can reasonably be extrapolated to suffi- 
ciently high-T and to Brillouin-scattering frequencies. In 
Section III this formalism is applied to available high 
quality sonic and ultrasonic data on u-SiOa, covering 
more than four orders of magnitude in the frequency v. 
The model parameters that are thus extracted are then 
used in Section IV to estimate the TAR contributions 
to both velocity dispersion and damping at Brillouin- 
scattering frequencies. The anharmonic contribution can 
then be extracted from the total signal. It is analyzed 
in terms of a mean thermal mode relaxation time r t h in 
Section V. Interestingly, this time is found to be about 
an order of magnitude longer than that of densified sil- 
ica glass, d-Si02- Section VI presents a synthesis of the 
above analysis for the entire range of v, from ultrasonic to 
Brillouin-scattering frequencies. A part of that analysis 
concerns the velocity changes with T. We find an anoma- 
lous increase of the bare velocity i>oo with T, nearly the 
same in t>-Si02 and d-Si02- This unrelaxed velocity, Voc, 
is obtained after subtraction of the velocity changes pro- 
duced by both thermally activated relaxation and anhar- 
monicity. The T-dependence of Woo suggests that silica 
experiences a progressive structural change with increas- 
ing T. Section VII is a discussion, which mentions the 
crossovers between the various frequency regimes, consid- 
ers the suitability of power laws for the description of the 
dependence of the damping on v, and evoques possible 
extension to other glasses. 



II. THE PHENOMENOLOGY OF THERMALLY 
ACTIVATED RELAXATION 

We consider an assembly of defects represented by 
double-well potentials that are separated by barriers of 
height V and whose depths differ by the asymmetry A. 
The energies V and A are randomly distributed accord- 
ing to a distribution P(A, V) dA dV to be discussed be- 
low. The system is thought to hop continuously between 
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the wells. The energy difference between the wells is cou- 
pled to the strain e of a sound wave of angular frequency 
f2 = 2nv by a deformation potential 7 = ±<9A/<9e. Ow- 
ing to the delayed energy exchanged in hopping, this pro- 



duces the dissipation of the sound wave. This situation 
has been described in great details elsewhere |MI(l lTllll2|| . 
To an excellent approximation, it leads to a relaxational 
contribution to the internal friction given by 



2 i>OC />OC A ( \ 

Q-A = J jA I dVP(A,V)s^- TT ^ I . (la) 

In this expression, p is the material density, v is the velocity of sound, and T is in energy units. We remark that both 
7 and v depend on the polarization of the acoustic wave, longitudinal (LA) or transverse (TA). In this writing, V is 
restricted to positive values while the distribution is symmetric in A. We also note that P(A, V) dA dV is a density, 
i.e. a number per unit volume. The associated velocity change, which follows from the Kramers-Kronig relation, is 
given by 
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In these equations, r is the relaxation time for hopping 
within the double well. It is given by 

V , A 

t = t exp — sech — , (2) 

where To is the inverse of an attempt frequency, as shown 
in detail e.g. in [l2T ]. 

The key in applying these expressions is to use a rea- 
sonable distribution P(A, V). For small A and V the 
distribution is often replaced by a constant P. This is 
suggested by the T-dependence of the specific heat at 
low T which only probes low values of A and V 0] • Of 
course, a constant P cannot be extended to high values 
of A and V as this leads to a diverging integral den- 
sity of defects which is unphysical. A reasonable guess 
for P(A, V) can be obtained with the help of the soft- 
potential model (SPM) [l4|. That model is character- 
ized by a distribution of random dimensionless cubic and 
quadratic coefficients, £ and 77 respectively, by an energy 
scale of the potential and by a characteristic cross- 
over energy W <C £q We are only interested here 
in the region 77 < with 1 ^> \rj\ ^> £ 2 which gives 
double wells with barriers centered at the origin of the 
soft mode coordinate x. As shown in ^5|, owing to the 
latter choice, the variables £ and 77 are not statistically 
independent. This leads to a sea-gull singularity in their 
distribution, P(£,rj) — ^\ri\Vo(^, rj), where Po is fi- 
nite near the origin. For the range of values of inter- 
est here, one has A = £ -i=£|77| 3 / 2 and V — ±£ M 2 - 
The deformation potential of the SPM is also function 
of 1 77 1 with 7ap M oc 1 77 1 as defined in and further ex- 
plained in [13. The terms 7 2 P(A, V) dAdV of (la) are 
transformed into 7spm-P(£> v) d£ dr] in the SPM. Using 
the Jacobian \d(A,V)/d(£,rj)\ = £q^\t]\ 5 / 2 , one finds 

that P(A,V) tx V-V^Voi^v) 0- This wa s already 
used by Keil et al. 0] who selected a distribution P(V) 



proportional to times a gaussian. These authors ex- 
perimentally found that indeed £ is very near 1/4 in the 
case of silica. 

For convenience, and for lack of different compelling 
indications, we assume that the distribution can be fac- 
tored into 

P(A,V) = f(A)g(V) . (3) 

For g(V), inspired by 0] and [l8|. we use the normalized 
form 

Since £ < 1, this expression is integrable. The norm N s 
is selected so that g(V) integrates to 1. This gives 

/•OO 

N- 1 = / x~ c exp(-^x 2 )dx = r(l-C) U(i-C,0) . 
Jo 

(4&) 

The function £/(a,0) is the parabolic cylinder function 
|19| |. As explained in the following Section, our indepen- 
dent analysis of a collection of data larger than that used 
in 0| also leads to £ very close to 1/4. This provides a 
solid support for this particular choice of g(V). For /(A) 
we use a simple gaussian, rather similar to the gaussian 
cut-off of the linear asymmetry coefficient D\ used by 
Gil et al. 20]. Indeed, D 1 cx £\r)\ oc AV' 1 ^ 4 , so that a 
gaussian in D\ is very close to one in A, the power of 
V connecting the two variables being quite small. This 
functional form for /(A) was already employed succes- 
fully by Bonnet [2l|. It is written 

/ 1 A 2 \ 

/(A) = 7V f / exp(~-^-J , (5a) 

where fo is defined by the normalization condition that 
the integral of /(A) equals /oVq- The norm JVf is then 
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dimensionless and given by 

/oo 
exp(-±5 2 x 2 )dx = V2n/S , (56) 
- OO 

where we defined 

S = V Q /Ac , (6) 

Ac being the cut-off value of the asymmetry. We finally 
remark that with the above definitions, the integral den- 
sity of defects is 

/OO />OC 
dA / dVP(A,V) = foV . (7) 
- oo JO 

I 



To complete the calculation, we now introduce (4) and 
(5) into (1), and use (2). In performing the integrals 
in (1), we make the same approximation as in |(| that 
sech(A/2T) is replaced by 1 for |A| < 2T and by zero 
otherwise. This eliminates the sech factors in (1) and 
(2) and simply replaces the limits of integration on A by 
±2T. It is convenient to define a dimensionless constant 



C = ^hNjpv 2 . (8) 



One obtains 
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where $(z) is the error function, 
2 f z 

$(z) = — / exp(-x 2 )dx . (10) 
V 71 " Jo 

For the purpose of comparison with literature results, it 
will be useful to relate the value of C with the "tunneling 
strength" defined by C = P-f 2 / pv 2 in the usual tunneling 
model [III Ho . l22| . Given the distribution (4), one cannot 
strictly define P, but this does not prevent defining C in 
a consistent manner [l^. We calculate this for £ = J. 

We use 7 | PM = ±A 2 \r]\ and C = A 2 T r/ 7 L /2 / Wpv 2 from 
[l6||. where 7/£ = W/£q. Introducing these definitions 
into 7spm p (C, if) d£ dr] = 7 2 P( A, V) dA dV, and using 
(8), we find in the limit of small A and V, 

[2 W 3 / 4 Vn /4 

C = J- - . C ■ (11) 

V 7T A c 

Information on the numerical handling of Eqs. (9) and 
our choice of suitable fitting parameters are found in Ap- 
pendix A. 

III. ANALYSIS OF SONIC AND ULTRASONIC 
RELAXATIONAL DATA 

It is generally agreed that sonic or ultrasonic damping 
at temperatures above ~ 10 K is dominated by ther- 
mally activated relaxation. To obtain the model param- 
eters Vq, Ac, and tq entering Eqs. (9), it is necessary 



to analyze acoustic results over a sufficiently large range 
of frequencies v, this up to high ultrasonic frequencies. 
It implies comparing data from various sources and gen- 
erally acquired with different measurement techniques. 
TAR leads to a peak in Q~ e \ as observed by Anderson 
and Bommel in ultrasonic pulse-echo measurements |l|. 
These authors report precise data obtained on the TA- 
mode of vitreous silica at v = 20 MHz. This is one of 
the curves that will be used for our analysis, as displayed 
in Fig. 1. Other high-frequency data were measured us- 
ing the Bragg diffraction of light, 23] in a set-up similar 
to the original Debye-Sears experiment |2^|. Damping 
results at v ~ 200 MHz on both LA- and TA-modes 
have been obtained in that manner |23j |. From these, 
and from results presented in [25j, it is clear that LA- 
and TA-modes lead to identical Q" 1 peak shapes. This 
implies that similar distributions of defects are active in 
the damping of all acoustic waves, independently from 
their polarizations. We use the data on the LA-mode 
at 207 MHz from [23|, also shown in Fig. 1. The inde- 
pendence from polarization allows including in the eval- 
uation sonic frequency results acquired on macroscopic 
vibrational modes, such as in vibrating reed measure- 
ments. We use the data at v — 11.4 kHz of Classen, as 
reported in |T^. This curve is also displayed in Fig. 1. 
Results at intermediate frequencies were collected using 
composite oscillators, at 660 kHz on the LA-mode in 
and at 180 kHz on a torsional mode in As remarked 
in |18|. this particular method can easily lead to instru- 
mental background- loss contributions. This might have 
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FIG. 1: The data points are internal friction results on u-SiCh 
taken from the literature, [23|. 0, and [T^l in order of 
decreasing frequency. The solid lines show our adjustment to 
Eq. (9a) as explained below. 



been the case in Q, as suggested from data on other 
glasses presented by the same authors in 26] which show 
long absorption tails at high T. This notion is also sup- 
ported by fits that are explained below. Hence, it is the 
180 kHz data from [3 which is included in Fig. 1. With 
these four curves, the analysed data covers more than 
four decades in v with a nearly linear progression in log v. 

The curves shown in Fig. 1 have been simultaneously 
adjusted to Eq. (9a), allowing for each measurement an 
independent coupling parameter C. In the Marquardt- 
Levenberg routine used to this effect, the weights have 
been adjusted so that each curve had approximately the 
same weight, independently from the number of data 
points at each particular frequency. The excellent re- 
sults of this fit, with C fixed at 1/4, are illustrated 
by the solid lines in Fig. 1. The distribution param- 
eters corresponding to these lines are Vo = 659 ± 19 
K, log 10 T = —12.2 ± 0.14 with t in seconds, and 
Vo/Ac = 8.2 ±0.6. If the same fit is repeated with £ free, 
one finds C = 0.28 ±0.03, and similar values for the other 
parameters: V = 667±21 K, log 10 r = -12.2±0.18, and 
Vo/Ac = 7.7±0.7. This shows that the distribution g(V) 
suggested by the SPM is very adequate indeed. It also 
shows that C can be fixed at 1/4, which will be done for 
the rest of this paper. The rather high value obtained 
for Vb/Ac emphasizes that one should not neglect the 
cut-off in the asymmetry distribution. We return to this 
point below. 

Our analysis reveals a certain difficulty in comparing 
the absolute size of the various curves. This can be re- 
lated to some extent to the different polarizations of the 
waves in these measurements, or also to a small part 
to the different qualities of the silica samples employed. 
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FIG. 2: The data points of the two lower curves are relative 
velocity changes of «-Si02 taken from the literature, for 660 
kHz and |2^| for 6.8 MHz. The solid lines show our calculation 
of the TAR contribution according to Eq. (9b) with C — 1.4 x 
10 -3 . The data points of the upper two curves show the 
experiment minus the TAR contribution. 



However, it seems more likely that it arises to a large 
part from calibration inaccuracies in some of these mea- 
surements. Indeed, discrepancies in the size of Q~ l are 
directly observed in comparing literature reports, for ex- 
ample for the TA- waves at 20 MHz in 1] and 25]. The 
difficulty is also seen in comparing the LA to TA values 



of Q- 



23| and 



While in the former case it is 



the attenuations that have the same amplitudes for both 
LA and TA, in the latter it is the internal frictions which 
are practically the same for both polarizations. This can- 
not be reconciled as attenuation and friction differ by a 
factor proportional to the sound velocity, and the latter 
is ~ 5900 m/s for LA and ~ 3800 m/s for TA- waves at 
low T. We also observed that the results at 660 kHz in 
are fitted extremely well with the same model param- 
eters, provided one allows for a relatively small constant 
background contribution to Q~ x , equal to 0.1 x 10~ 3 . 
The values of C obtained for these various measurements 
range then from 1.1 x 10~ 3 in the case of [|| to 1.8 x 10~ 3 
for [J. These variations in C do not seem correlated to 
either the measuring frequency or the wave polarization, 
which is another reason to suspect calibration difficulties. 

It is necessary to obtain a more reliable value of C 
to analyze the Brillouin scattering data in the following 
Section. To this effect we remark that velocity measure- 
ments are generally both accurate and precise, allowing 
to follow small changes in Sv/v with excellent accuracy. 
We use two measurements of the LA-wave velocity, one 
at 660 kHz from I™ and the other at 6.8 MHz from 



|23| . The experimental values (Sv/v) eKp are shown in 
Fig. 2, together with continuous lines marked (Sv/v) Ie i. 
These are calculated from the relaxation equation (96) 
using the above values of Vq, Ac, and To, together with 
C = 1.4 x 1CP 3 . With the latter value, the dip around 
50 K is completely accounted for. This is emphasized by 
the difference (5v/v) exp — (6v/v) rc \ also shown in Fig. 2. 
Hence, we adopt for the following the coupling constant 
C = 1.4 x 10~ 3 which happens to fall within the range 
of values obtained by fitting the various Q _1 curves dis- 
cussed above. 

We now return to the small cut-off value Ac — 80 K 
which we obtained. This aspect has mostly been ignored 
by other workers. It was shown in j2l| that it is nec- 
essary to include a cut-off to obtain a good fit to the 
high-T tails in Fig. 1. For example, forcing Ac = Vq, the 
quality of the fits degrades appreciably above ~ 70 K, 
especially for the highest v curve. More importantly, the 
calculated values of (5v/v) re i are then much too large to 
account properly for the dip around 60 K observed in 
Fig. 2. Obviously, by the virtue of the Kramers-Kronig 
transform, (5v/v) Te \ integrates over a large part of the 
distribution, which is the reason for this problem. One 
might wish to gain an intuitive picture of why Ac can 
be so much smaller than Vq. To this effect, one can con- 
sider the models of TAR drawn in Fig. 2 of 27] . We 
take as simplest examples model A, in which the two Si 
atoms of a Si— O— Si bond are too close, and model B, 
in which they are too far apart. In either double 
well potential for the connecting oxygen results with a 
barrier height V that comes mainly from the separation 
of the two Si. This barrier can thus be quite high. On 
the other hand, the asymmetry is produced by the differ- 
ence in the wider environment of the two wells. In a hard 
glass, one might expect that these environments, which 
are dictated by the minimization of the energy, can be 
mostly quite similar. From these considerations, one in- 
tuitively anticipates that the ratio Vb/Ac might depend 
significantly on the particular glass. 

Finally, we remark that the temperature T max of the 
peak positions in Fig. 1 depends linearly on lnf2. A sim- 
ilar observation was already made in |27| . This typi- 
cal Arrhenius behavior supports TAR as the principal 
relaxation mechanism to describe Fig. 1, as opposed 
for example to incoherent tunneling that would lead to 
T max tx Vti as seen from Eq. f2.95) of pf. 



IV. ANALYSIS OF BRILLOUIN SCATTERING 
RESULTS 

High-resolution Brillouin scattering measurements of 
the temperature dependence of the LA-linewidth of vit- 
reous silica have been reported in [2jJ and 0] . Both ex- 
periments were performed near backscattering and below 
room T, using as exciting radiation the blue argon-laser 
line at Al = 488 nm. Measurements above room T are 
reported in |30j . These were performed at Al = 514 nm. 
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FIG. 3: a: The data points are Brillouin shifts obtained on 
u-SiCh near backscattering, from pfflj and pffl^ . The values of 
30] have been rescaled to the blue argon-laser line, Al = 488 
nm. The solid line is a guide to the eye. b: The internal 
friction derived from the Brillouin linewidths. The crosses 
are from |29|| and [3C|. while the dots are from [T^|. The solid 
curve shows the TAR contribution calculated from Eq. (9a) 
using the value of C determined from the velocities in Fig. 2. 



The Brillouin frequency shifts vb from [2^] and [3(j, the 
latter rescaled to Al = 488 nm, are displayed in Fig. (3a) 
over the entire range of T. 

The frequency widths extracted from the Brillouin 
peaks in |29| and |30| , after correction for the instrumen- 
tal broadening, have been converted to internal friction 
as shown in Fig. (36). The internal friction reported by 
Tielburger et al. in [12j is also shown on the same figure. 
One notices the remarkable agreement between these two 
independent determinations. It emphasizes that Bril- 
louin scattering gives an absolute measurement of Q , 
independent from calibrations, since it suffices to take 
the ratio of the Brillouin full-width to the Brillouin shift 
to extract Q -1 . This statement of course assumes that 
the spectroscopy can be performed with the required res- 
olution, and in particular that the broadening produced 
by the finite collection aperture can be made sufficiently 
small. 

Also shown in Fig. (36) is the curve Q~ el calculated 
with Eq. (9a) for v = 35 GHz and with the parameters 
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FIG. 4: The lower data points show the measured ultrasonic 
velocities from Fig. 2 together with the hypersonic velocity 
derived from Fig. (3a). The upper data points show the same 
velocitites after subtraction of the TAR contribution. We note 
that after this subtraction there remains a small difference 
between the slopes observed at high T for the two different 
frequencies. These slopes are shown by the dashed lines and 
they are discussed in Appendix B. 



determined in Section III. It is evident that at Brillouin 
frequencies TAR cannot account for the entire internal 
friction. The value C ~ 1.4 x 10~ 3 is confirmed by the 
Brillouin velocity results explained below. However, even 
if one would allow oneself to arbitrarily double the cou- 
pling constant C, the line in Fig. 3(6) would still not 
superpose the data points. In particular the peak in Q~J 
occurs at a lower T- value than the peak in the observed 
Q -1 . The difference between the two curves is attributed 
to its largest part to the anharmonic coupling of the LA- 
waves with the thermally excited vibrational modes, as 
reported in (I) for densified silica glass. 

The velocities are extracted from the Brillouin frequen- 
cies i/b of Fig. (3a) using v = Al^b /2nsin|, where n is 
the refractive index and 9 is the internal scattering angle, 
here close to 180°. To derive precise values on v(T), it is 
necessary to know n(T) with an equal precision. This in- 
formation has been derived from 31] . The results below 
300 K are shown in Fig. 4, together with the velocities 
at 6.8 MHz from [2j|. The relative changes reported 
in [|| have also been converted back to velocities using 
v = 5888 m/s. Subtracting from all three curves the 
(6v) re \ calculated from Eq. (96), one observes that be- 
low 150 K all the data collapse quite well. This confirms 
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FIG. 5: The internal friction observed in Brillouin scattering 
after subtraction of the TAR contribution shown in Fig . (36). 
The crosses are from |2fl| while the dots are from |l2l |. The 
solid line is a rough estimate for the small tunneling contri- 
bution extrapolated from data presented in |3^l . 



that the value of C is also correct at Brillouin frequen- 
cies. It also emphasizes that the ^-dependence in the 
depth of the dip around 70 K is well predicted by TAR. 
As explained in Section VIB, the velocity changes that 
are produced by the anharmonicity, (<5u) a nh, show little 
dispersion compared to the large dispersion in (Sv) re \. 
Hence, they do not modify the above conclusion. 

Figure 5 shows the difference between the experimen- 
tal internal friction and the calculated TAR contribution. 
Similarly to the results on d-Si02 in (I), there is a region 
where the signal falls rapidly with decreasing T. How- 
ever, it does not seem to fall to zero sufficiently fast. 
In spite of the large scatter in the experimental values, 
there is a hint for an additional contribution to Q^ 1 at 
low T, in the region from ~ 20 to ~ 60 K. At these 
hypersonic frequencies and temperatures, one can antic- 
ipate a contribution Q^ls arising from the relaxation of 
two-level systems (TLS) by incoherent tunneling. Gen- 
erally, Qtls(^) exhibits a plateau separating a regime 

^Tmin ^ 1 at low T <T from fh" m i n < 1 at higher T, 

as shown in [28| where r m ; n and Toe fi 1 / 3 are defined. 
The position, extent, and height of this plateau are in- 
dependent, as emphasized in Fig. 6 of (32]]. According 
to |23, the extent of the plateau should shrink to zero 

at sufficiently high frequencies such that T> T*, where 
T* is a frequency independent characteristic tempera- 
ture marking the end of the plateau region. In w-Si02, it 
has been observed that T* ~ 5 K 0. The plateau be- 
comes then a hump slightly above T- As shown in [32| . 
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for v — 46 MHz one should have T— 3 K. This value 
increases with fi 1 / 3 , as confirmed from sound-damping 
data between 330 and 930 MHz |3^, showing a shoul- 
der on the low T side of the TAR peak. Hence, at 35 

GHz one expects 30 K, giving a broad TLS hump 
centered around 40 K. Q^ s oc T 3 below this hump, and 
Qtls x r ^~ 1 a bove it |2g. Extrapolating the observa- 
tions reported in we posit 

_ 1.3 x 10^3 
^ TLS ~ (30/T) 3 + (T/30) 1 > 

as a rough estimate for this relatively small contribution 
to Q^ 1 . This curve, which peaks around 40 K, is shown 
in Fig. 5. The difference 

Qanh = Q7 e ] QtL (13) 

will be used below to analyse the anharmonic damping. 
Since there seems to be a small tunneling contribution 
Qtls ' then there ought to be by the Kramers-Kronig re- 
lation a corresponding contribution (6v)tls- As shown 
in Eq. (2.96) of j^, this should produce at high-T a 
frequency dependent term in Tin CI. As shown in Ap- 
pendix B, there are reasons to believe that it is this term 
that produces the difference in slopes indicated by short 
dashed lines in Fig. 4. 

To conclude, we find a relatively large a t hyper- 
sonic frequencies. It is of the same order of magnitude 
as Q~ c \- The small additional term Q^ls th a t nas been 
discussed above is a correction of minor importance com- 
pared to Q7 n hi as obvious from Fig. 5. 

V. THE THERMAL RELAXATION TIME 

In the spirit of (I) we now analyze the anharmonic 
damping with the expression 

«i-^rra^- < 14 » 

where r t h is the mean lifetime of the thermal modes and 
A(T) = j 2 C v Tv/2pvl is the prefactor given in Eq. (56) 
of (I), with C v the specific heat per unit volume and 
i>d the Debye velocity. We remark that all quantities 
entering A(T) are known, except for the mean-square 
average Gruneisen parameter r ) 2 . Thus we can directly 
plot the quantity 1 2 Q^^/A as it is independent from 7 2 . 
The result obtained using the values taken from 

Fig. 5 is drawn in Fig. (6a). Points below 60 K are not 
shown as the data is too noisy at low T to be significant 
after division by A cx C V T. This plot exhibits a peak 
for firth = 1- From the height of this peak one reads 
7 2 ~ 3.6. With this determination, Eq. (14) can then 
be solved for firth, a nd thus for r t h- The procedure is 
explained in Appendix C. The result is shown in Fig. (6b). 
One observes that 7th (T) is proportional to T -2 over the 
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FIG. 6: a: The anharmonic damping contribution scaled by 
A(T)/^ 2 . The solid line is a guide to the eye calculated from 
the corresponding line in b. The peak is located at Tmoi (see 
Appendix C) and its height gives 7 2 /2. b: The values of r t h 
for w-Si02 obtained from the data in a after solving Eq. (14). 
Data for x-cut crystal quartz and d-Si02 from (I) are shown 
for comparison. The lines for «-Si02 and d-Si02 are guides 
to the eye which are sums of negative powers of T, while the 
line for c-Si02 is calculated from Eq. (15). 



main part of the data, from ~ 100 K to ~ 300 K. At 
higher T, r t h tapers off, while at lower T there is a hint 
that Tth(T) might diverge faster than in T~ 2 . The solid 
line is a guide to the eye which is a sum of powers of T. 
The same law is then used to trace the guide to the eye 
in Fig. (6a). 

In Fig. (66) we have traced for comparison the values 
of the relaxation time of the thermal modes for two other 
cases. Firstly, we observe that 7th (T) obtained in (I) for 
densified silica glass, d-Si02, is substantially shorter than 
in w-Si02- However, it shows a very similar functional 
dependence on T. Secondly, in order to deepen the un- 
derstanding, we applied the same method of analysis to 
data on the LA-mode of X-cut crystal quartz, c-Si02- Wc 
used to this effect the attenuation coefficient measured by 
Bommel and Dransfeld at 1 GHz [34[. At temperatures 
above 150 K, we completed these data using the width 
of the Brillouin line reported in |3^. These data were 



analyzed similarly to Fig. (6a) leading to j 2 ~ 0.95. The 
corresponding T t h(T) is shown in Fig. (66). In this case, 
the line through these points is not a guide to the eye as 
for the two glasses, but it is an independent determina- 
tion of 7th(r) using the well-known kinetic expression for 
the thermal conductivity, 

K = | CvWoTth • (15) 

We observe that in this case the value of rth derived 
from k (line) is in remarkable agreement with that ob- 
tained from the measurement of anharmonic damping of 
the acoustic modes. This emphasizes the physical sig- 
nificance of T t h- it really is the mean lifetime of the ex- 
citations in the thermal bath which at sufficiently high 
T is entirelly controlled by Umklapp (U) processes. As 
pointed out in f36:j, this produces a relaxation propor- 
tional to the phonon population, and thus at high T one 

has r t h oc T _1 . At intermediate T, where C v oc T and 
k cx T _1 , Eq. (15) indicates that r t h oc T~ 2 . This is 
indeed observed in c-Si02 over a large range of T, as seen 
in Fig. (66). 

Equation (15) assumes propagating thermal phonons. 
For this reason, it becomes invalid in glasses at soon as 
T increases beyond the thermal conductivity plateau lo- 
cated around 10 K in t>-Si02. However, we remark that 
the functional dependence Tth (J 1 ) observed in the two 
glasses is very similar in shape to that in c-SiC>2. We 
have no simple explanation for this. In glasses, there are 
two competing effects that modify the picture presented 
above. On the one hand, the strict quasi-momentum 
conservation which is invoqued in U-processes is strongly 
relaxed, and this must greatly enhance the interactions 
in the thermal bath, decreasing Tth- This could account 
for the observed difference between crystal and glasses. 
On the other hand, the thermal modes are not expected 
to be propagating plane-waves but they are at best dif- 
fusive. This restricts the spatial extent of the modes, 
greatly decreasing their overlap and thereby their inter- 
actions, which increases T t h- We believe it is the lat- 
ter effect which produces the much longer Tth i n v-Si02 
compared to d-Si02- Indeed, in d-Si02 the boson peak 
is strongly reduced [33 ■ This increases the crossover fre- 
quency lo co beyond which the acoustic excitations become 
diffusive, as recently confirmed by inelastic x-ray scatter- 
ing [38j. This in turn increases the mean spatial extent 
of the modes and thereby their interactions leading to a 
faster thermalization. 

Finally, we remark that the above analysis gives val- 
ues for the mean Gruneisen parameter 7 2 . We find 0.95 
for the LA- mode of x-cut quartz, 3.6 for w-Si02, and 8 
for d-Si02- In the case of crystal quartz, the agreement 
between the line calculated from the thermal conductiv- 
ity k and that obtained from Q _1 gives a solid support 
for the value of 7 2 . Unfortunately, one cannot make a 
similar comparison for the glasses. 
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FIG. 7: The internal friction at 35 GHz from the Brillouin 
scattering measurements of [13 (o) and H| (+). The solid 
curve is the entire adjustment made here which is the sum of 
three contributions Q~f, QtlS' an< ^ Qlnh- 



VI. SUMMING-UP 
A. The internal friction 

Figure 7 summarizes the present analysis of the in- 
ternal friction observed in U-S1O2 with Brillouin scatter- 
ing around 35 GHz. One identifies two main contribu- 
tions: Q~J arising from thermally activated relaxation 
and Q~nh = Q^ 1 - (Qrel + Qtls) which results from 
network viscosity. Around room T, Q^nh ^ s nearly twice 
as large as Q~J . The incoherent tunneling contribution 
*5tls ' h° h really exists at these high frequency and tem- 
peratures, mainly produces the small hump around 40 K. 
Although there is a hint for such a feature in the data 
shown in Fig. 7, it would need to be confirmed by more 
precise measurements. Also, we neglected in this analysis 
a possible contribution arising from the quasi-harmonic 
oscillators of the soft-potential model This is ex- 

plained in Appendix D. 

As v is increased beyond 35 GHz, we find that the value 
<5~J calculated with Eq. (9a) saturates. This saturation 
essentially results from the cut-off in /(A). For example 
at room T, Q~J does not increase beyond ~1.5xl0~ 3 . 
On the other hand, at sufficiently high T, Q~i contin- 
ues to grow with f2 according to Q~h = ^^ r th which 
applies as long as SlT t h <C 1. We also find that Arth is 
practically constant from 100 to 300 K. For these reasons, 
the total Q _1 becomes dominated by network viscosity, 
it continues to increase proportionally to Q up to several 
hundred GHz, and it is nearly T-independent from 100 K 
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FIG. 8: a: Calculated velocity changes produced by TAR (rel) 
and by network viscosity (anh) at two typical frequencies, b: 
The measured velocities at the same frequencies (lower data 
points) and the unrelaxed values v x after correction for TAR 
and network viscosity (upper data points). The line is a fit of 
Voo to Wachtman's equation as a guide to the eye. The Voo 
data for d-Si02 from (I) are also shown for comparison. 



to 300 K. This is in excellent agreement with the results 
on acoustic propagation reported for u-Si02 thin films by 
Zhu, Maris, and Tauc j^. These authors find that the 
mean free path, £ _1 , grows in Jl 2 from —30 GHz to —300 
GHz, and that it is nearly independent of T from 80 to 
300 K. Since l~ x = Q^Sl/v, the dominance of Q^nh> ^ ts 
growth oc Q, and its near constancy in T, fully account 
for the data reported in (3^. Our predictions are also 
fully consistent with three data points measured with 
UV Brillouin scattering from -50 to -100 GHz pol lllf . 

On the other hand, at ultrasonic frequencies, i.e. much 
below 35 GHz, the contribution of which decreases 
proportionally to SI becomes completely negligible com- 
pared to Q~j- This justifies the analysis based only on 
TAR that was performed in relation with Fig. 1. 



B. The velocity changes 

Figure (8a) shows the velocity changes calculated from 
Eq. (96) for TAR, and using 

f 6v \ A 1 

l»L a -^ (16) 

for the network viscosity, this at two typical frequen- 
cies. One notices that the increase in v from 6.8 MHz 
to 35 GHz leads to a large reduction of |<5v rc i|. On the 
other hand, |5w an h|, although quite large above 100 K, 
changes relatively little with v. In particular, the dashed 
curve is entirely in the regime firth <C 1, so that it is prac- 
tically independent of SI at lower Sl-values. This justifies 
the approach used to determine C in Fig. 2 where the an- 
harmonicity contribution was simply neglected. At suf- 
ficiently high-T, anharmonicity always leads to a quasi- 
linear decrease of v with T, very similar to what is ob- 
served in crystals 0> 03 • This has been seen in ultra- 
sonic measurements in a large number of glasses, e.g. in 
[44) 1 where this effect is interpreted as such. At these fre- 
quencies, one is in a regime where the velocity changes 
are dominated by network viscosity while the damping 
mainly originates from TAR. It is generally necessary to 
go up to Brillouin frequencies to observe in Q" 1 both 
contributions simultaneously. 

The velocity changes in tetrahedrally coordinated 
glasses, in particular in v-Si02 and in d-Si02, are more 
complex than described above. Fig. (86) shows the ob- 
served velocities, v, and the unrelaxed value that is cor- 
rected for both TAR and network- viscosity dispersions, 
Voo = v — [(Sv) Te i + (5w) a nh]- As observed, the exper- 
imental points collapse to a which on this scale is 
essentially Sl-independent. However, rather than being 
constant, increases considerably with T. In the same 
figure, the results obtained in (I) for d-Si02 are shown 
for comparison. Both solid lines are guides to the eye 
obtained by an ad hoc adjustment to Wachtman's equa- 
tion |45| . The behaviour of is remarkably similar in 
both glasses, pointing to an anomalous hardening of silica 
with increasing T. As discussed in (I), the earlier pro- 
posal that this is a manifestation of large structural in- 
homogeneities [4.6J] is not supported by observations. The 
anomalous hardening has recently been attributed to a 
progressive local polyamorphic transition associated with 
abrupt rotations of randomly distributed Si-O-Si bonds, 
found in simulations |47| . Thermal agitation would redis- 
tribute the Si-O-Si bond angles in configurations where 
they are more resilient. However, contrary to |4S| . our re- 
sults indicate that a very similar mechanism would then 
be active in permanently densified silica. Repeating the 
warning made in (I), although has a clear physical 
significance, it would be nearly impossible to perform at 
the moment an experiment to directly measure it. This 
is even more so in u-Si02 than in d-Si02, owing to the 
lower uj co of the former glass. 
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C. The tunneling strength 

It remains to compare the value C = 1.4 x 10~ 3 found 
in Section III to the tunneling-strength parameter C ~ 
3.0 x 10 -4 from the literature. To do this, both values 
are inserted in Eq. (11) to calculate the cut-off energy 
W that applies to the present case. We find W ~ 7 K, 
as opposed to the TLS value W ~ 4 K derived from the 
low-T specific heat data 0|. Such a larger value of W 
was already reported, resulting from a similar analysis 
in [lfl| . Its origin was convincingly explained there as 
arising from a different weighing of the distribution of 
W- values, in W~ A for the specific heat, and in W^ 1 for 
the ultrasonic absorption. The value of C is thus fully in 
line with expectations. 

VII. DISCUSSION 

We have shown that ultrasonic and hypersonic damp- 
ing and dispersion in vitreous silica at temperatures 
above the quantum regime originate in at least two main 
processes: thermally activated relaxation and network 
viscosity. These have different f2 and T dependences 
which leads to a non-trivial variation of the scattering 
linewidth T with f2. At a given T, TAR dominates the 
damping at low Q, while anharmonicity dominates at suf- 
ficiently high fl. Between these two regimes there occurs 
a "dynamical crossover" which for v-Si02 at room T falls 
near our Brillouin frequency, as shown in Fig. 10 of |39| 
and confirmed in Fig. 7 above. Above this crossover, the 
damping is dominated by Q^nh = ^^ r th; which implies 
that r oc ft 2 . As ft is increased far beyond optical Bril- 
louin frequencies, another broadening mechanism eventu- 
ally comes into play. Around 300 GHz, one expects the 
onset of a new contribution to Q~ x , growing with a higher 
power of ft, possibly in ft 3 , as suggested in This ad- 
ditional contribution is currently thought to result from 
the hybridization of the acoustic modes with boson-peak 
excitations |^fj, |&lj , as inelastic x-ray scattering observa- 
tions on c?-SiC>2 strongly suggest |33,|53. Unfortunately, 
inelastic x-ray scattering at sufficiently small scattering 
vectors and with sufficiently high energy resolution could 
not yet be performed on v-Si02 to settle the issue of this 
onset. This hybridization eventually leads to the Ioffe- 
Regel crossover at uj co , above which the excitations have 
lost their plane-wave character. Then, the wave-vector q 
does not remain a good label for the modes which are at 
best diffusive |5^|. The excitations of the thermal bath 
are expected to be of that nature in glasses, as discussed 
in Section V. 

Power laws T oc fl n , or T oc q n which is equivalent as 
long as the frequency is well below the Ioffe-Regel limit 
so that q oc ft, are often employed to represent the damp- 
ing data r(fi). As just explained, r oc ft 2 is a reason- 
able approximation in the region between the dynami- 
cal crossover and the Ioffe-Regel one. On the opposite, 
such laws should be viewed as quite rough approxima- 



tions in the TAR-dominated regime. The reason is that 
the dependence of T on ft predicted by Eq. 9a is more 
complicated. It cannot be written simply in terms of a 
T-dependent exponent n(T). For example, in v-Si02 at 
room T, a value n = 1.3 is given in |54f. while n = 1.8 
is derived from [7j. We conclude that power laws in the 
TAR-dominated regime are ad hoc devices which only 
represent approximately the data, especially if the latter 
are sufficiently noisy to mask the disagreements. 

At this point, it is legitimate to ask to what extent 
anharmonicity might be important to sound damping in 
other glasses. The answer to this question depends on 
two factors: i) the strength, density, and distribution of 
the relaxing defects; ii) the strength of the anharmonic- 
ity. It is conceivable that in glasses that contain a large 
quantity of defects up to high values of Vq and Ac , TAR 
would mask the anharmonicity. In that case the dynam- 
ical crossover might move up to nearly uj co , so that a 
regime T oc ft 2 might not be observable at all. We have 
found such a case in lithium diborate, V-LI2B4O7 |5l| . 
This type of situation is expected for most polymers in 
which the tunneling strength C is quite large |22j . On 
the other hand, there exist many glasses or amorphous 
materials for which C is quite small, well below 10 -4 |22j| . 
Such are the amorphous group IV semiconductors, C, Si, 
and Ge. In that case, the anharmonicity is expected to 
dominate over a considerable range of ft and a calcula- 
tion to that effect has already been performed [E^. Of 
the tetrahedrally coordinated glasses, it will be partic- 
ularly interesting to reinvestigate u-Ge02 and u-BeF2. 
Difficulties met in obtaining a quantitative description 
of damping and dispersion in germania on the basis of 
TAR alone, described in [ljj], might find their resolu- 
tion by including in the analysis the cut-off Ac and the 
network-viscosity contribution. 

To summarize, TAR essentially explains acoustic 
damping results in silica glass up to a few GHz. At higher 
frequencies, the anharmonic coupling of hypersound to 
the modes that form the entire thermal bath becomes 
progressively dominant. This conclusion, in line with 
previous observations and simulations [55|, corrects 
statements to the contrary often found in the literature, 
e.g. in The anharmonicity which in our view 

produces Q~^ h 1S a coupling of sound to the whole bath, 
in addition to the relatively small fraction of modes form- 
ing the boson peak. This process, which we call network 
viscosity, is in a way quite similar to the usual lattice 
viscosity of crystals. 



APPENDIX A: ESTIMATING EQS. (9) 

For the Marquardt-Levenberg adjustment of Eqs. (9), 
it is convenient to use the parameters C, Vq, logioTo, C, 
and 5. Taking as integration variable x = V/Vq, and 
defining the auxiliary variables a = Vo/T and y = firo, 
Eq. (9a) becomes 
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Q-] = C${V26/a)a / a;- c exp(-i x 2 ) - 
Jo 1 



y exp(aa;) 
f y 2 exp(2ax) 



dx 



(Al) 



Good results are obtained with a simple trapezoidal in- 
tegration, taking a geometric progression for the vector 
of ir-values. We used x starting at ~ 10 -4 and ending at 
~ 10, in steps of 1%. Eq. (96) is handled similarly. 
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APPENDIX B: VELOCITY CHANGES DUE TO 
INCOHERENT TUNNELING 

One observes on the upper curves of Fig. 4 that there 
is a small difference in the high temperature slopes for 
the two different frequencies. The origin of this ef- 
fect might be in the velocity dispersion associated, by 
the Kramers-Kronig transform, with the relatively small 
damping Q^ls produced by incoherent tunneling. At 
high T, the damping is given by the third equation (2.94) 



The corresponding <5vtls is derived from the third equa- 
tion (2.96) of [2S|. which gives for the ^-dependent part 



Svtls - 0.45^ (lnfi)T 



(B2) 



The temperature T* marks the end of the low-T plateau 
region in Q^ 1 and it is given by Eqs. (2.78) and (2.79) 
of |28|. Owing to the logarithmic dependence in (B2), 
it is the difference of slopes at two frequencies which is 
significant, 



J 



A(5v T ls) _ (^i)tls - {5v 2 )tls 



T 



T 



0.45^In(fii/n 2 ) 



(B3) 



With the values in Fig. 4, and calculating T* following 
[28| . one finds C»7x 10~ 5 . As expected, this is smaller 
than the tunneling strength found from specific heat 
measurements at low temperatures [22j, the reduction 
being of the order of Aq/T. Hence, it is not unreason- 
able to assign the difference in slopes observed at high-T 
in Fig. 4 to the Kramers-Kronig transform of Q^ls- 



APPENDIX C: SOLVING EQ. (14) FOR THE 
RELAXATION TIME 

It is in principle trivial to solve the quadratic equation 
y = x/(l + x 2 ) for x = firth in terms of the measured 
y = Q^ 1 1 A. However, imaginary roots do occur in the 
region around x — 1 owing to fluctuations in the data 
leading to some points for which y > 1/2. These are the 
points whose ordinate lies above the maximum of the 
solid curve drawn in Fig. (6a). For these data points, 
we simply replaced x by 1, as the nearest real solution. 
This produces a line of data points with r t h = I/O in 
the presentation of Fig. (66). For the data points with 
y < 1/2, one must select between the upper and the 
lower root of the quadratic equation. This is set by the 
location of the maximum in the solid curve, Tjviax- For 
T < Tmsx-, the upper root applies since firth > 1, while 



the lower root applies at higher temperatures. 



APPENDIX D: ANOTHER POSSIBLE SOURCE 
OF DAMPING 

We have not considered the damping arising from 
the weak anharmonicity of the nearly harmonic oscil- 
lators (HO) of the soft potential model in the domain 
i] > This contribution should be most active 

at "intermediate" frequencies. In v-Si02, its strength 
has actually been pred icted to be largest near our 
Brillouin frequency |17| . It must be noted that the 
entire (O, T)-domain covered by the propagating pulse 
measurements of j3^| falls within this "intermediate" 
frequency range where the additional damping obeys 
Qho « CT/Vn [13 . If the strength of Q H q would be 
sufficient, one should definitely notice its effect on the 
(TL T)-dependence of the mean free path reported in 
|39j from 80 to 300 K and from 30 to 300 GHz. On 
the contrary, the mean free path is found practically 
T-independent and it increases oc f2 2 , in agreement 
with the network-viscosity contribution. It seems thus 
justified to neglect Qjjq m our analysis of the Brillouin 
results up to 300 K. The difficulty in estimating the 
size of Qho apparently lies in finding the proper value 
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for the tunneling-strength parameter C that applies 
to it. The value C = 3.0xl0 -4 , appropriate for TLS 
at very low temperatures [2^ . is certainly too large 
here. Indeed, using this value in Eq. (3.18) of [l^, we 
calculate Q^o - 5xlCT 5 x T(K), which at 300 K is 
larger than the entire Q observed. ^From an analysis 
of heat-release measurements there are grounds to adopt 
here a value of C which is at least five times smaller 1561 . 



Even so, the resul ting Qjjq is still too large to agree 
with the results of 39]. Our Brillouin measurements up 
to 700 K shown in Fig. (36) suggest that there could be 
a small contribution oc T that starts being felt above 
~ 400 K. This contribution would then be of the order 
of 10~ 6 x T(K), which implies a very small effective C 
for this particular damping mechanism. 
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